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Abstract 



We consider a variant of the standard electroweak theory in which the 
Higgs sector has been modified so that there is a classically stable weak scale 
soliton. We explore fermion number violating processes which involve soliton 
decay. A soliton can decay by tunnelling under the sphaleron barrier, or the 
decay can be collision induced if the energy is sufficient for the barrier to be 
traversed. We present a classical solution to the Minkowski space equations 
of motion in which a soliton is kicked over the barrier by an incoming pulse. 
This pulse corresponds to a quantum coherent state with mean number of 
W quanta ~ where g is the SU{2) gauge coupling constant. We also 

give a self-contained treatment of the relationship between classical solutions, 
including those in which solitons are destroyed, and tree- level quantum am- 
plitudes. Furthermore, we consider a limit in which we can reliably estimate 
the amplitude for soliton decay induced by collision with a single 1^-boson. 
This amplitude depends on g like exp(— c^"^/^), and is larger than that for 
spontaneous decay via tunnelling in the same limit. Finally we show that 
in soliton decays, light SU{2)l doublet fermions are anomalously produced. 
Thus wc have a calculation of a two body process with energy above the 
sphaleron barrier in which fermion number is violated. 
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I. INTRODUCTION 



In the standard electroweak theory, fermion number violation is present at the quantum 
level but these processes are seen only outside of ordinary perturbation theory. A baryon 
number three nucleus can decay into three leptons. The process is described as an instanton 
mediated tunnelling event leading to an amplitude which is suppressed by exp(— STr^/gf^), 
with g ^ 0.65 the SU{2) gauge coupling constant. At energies above the sphaleron barrier 
10, fermion number violating processes involving two particles in the initial state are gener- 
ally believed to be also exponentially suppressed (At energies comparable to but below 
the sphaleron barrier, Euclidean methods @] have been used to show that the exponential 
suppression is less acute than at lower energies, but the approximations used fail at energies 
of order the barrier height and above.) Unsuppressed fermion number violating processes are 
generally believed to have of order 47r/(yf^ particles in both the initial and final states. This 
all suggests that fermion number violation will remain unobservable at future accelerators 
no matter how high the energy, whereas in the high temperature environment of the early 
universe such processes did play a significant role 0. 

In this paper, we explore the robustness of these ideas by studying a variant of the 
standard model in which the amplitudes for certain fermion number violating collisions, as 
well as for spontaneous decays, can be reliably estimated for small coupling g. The model is 
the standard electroweak theory with the Higgs mass taken to infinity and with a Skyrme 
term ^ added to the Higgs sector. With these modifications, the Higgs sector supports 
a classically stable soliton which can be interpreted as a particle whose mass is of order 
the weak scale 0. Quantum mechanically, the soliton can decay via barrier penetration 
||8HlO|. Classically, i.e., evolving in Minkowski space using the Euler-Lagrange equations, 
the soliton can be kicked over the barrier if it is hit with an appropriate gauge field pulse. 
Correspondingly, the soliton can be induced to decay quantum mechanically if it absorbs 
the right gauge field quanta. Regardless of whether the decay is spontaneous or induced, 
ordinary baryon and lepton number are violated in the decay. We shall see that the model 
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has a limit in which fermion number violating amplitudes can be reliably estimated both 
for processes which occur by tunnelling and for those which occur in two particle collisions 
between a soliton and a single H^-boson with energy above the barrier. 



A. The Model 

To modify the standard model so that it supports solitons, proceed as follows. Note that 
in the absence of gauge couplings the Higgs sector can be written as a linear sigma model 

£h = ^ Tr -ji^ [*^*] - ^')' (1-1) 

where 

{ipo, ipi) is the Higgs doublet, and v = 246 GeV. One advantage of writing the Lagrangian 
in this form is that it makes the SU{2)l x SU{2)ji invariance of the Higgs sector manifest. 
The scalar field $ can also be written as 

$ = (7C/ (1.3) 

where U is SU{2) valued and cr is a real field. In terms of these variables 

1 / 2\ 2 

Cn^-a'Tr[d^U^d^u]+d^ad''a-xia'-^\ . (1.4) 



The Higgs boson mass is y2\v. We work in the limit where the Higgs mass is set to infinity 
and a is frozen at its vacuum expectation value v/\/2. Now 
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Ch = ylV [d^U^d^u] (1.5) 



which is the nonlinear sigma model with scale factor v. We will consider only those configu- 
rations for which the fields approach their vacuum values as |x| oo for all t. We can then 
impose the boundary condition 

hm t/(x,t) = l, (1.6) 

|x|— +00 



which means that at any fixed time U is a map from into SU{2). These maps are 
characterized by an integer valued winding number which is conserved as the U field evolves 
continuously. However if we take a localized winding number one configuration and let 
it evolve according to the classical equations of motion obtained from ( p..5| ) it will shrink 
to zero size. To prevent this we follow Skyrme ^ and add a four derivative term to the 
Lagrangian. The Skyrme term is the unique Lorentz invariant, SU{2)l x SU(2){i invariant 
term which leads to only second order time derivatives in the equations of motion and 
contributes positively to the energy. 

,2 _ _ 1 



32e2 



rTr 



[1.7) 



where e is a dimensionless constant. 

Of course this Lagrangian is just a scaled up version of the Skyrme Lagrangian which has 
been used P,|lT|,|12[ to treat baryons as stable solitons in the nonlinear sigma model theory 
of pions. To obtain the original Skyrme Lagrangian replace v in ( |1 . 7| ) by f^^ = 93 MeV. The 
stable soliton of this theory, the skyrmion, has a mass of TS/tt/c and has a size ~ 2/(e/7r) 
|T^. Best fits to a variety of hadron properties give e = 5.5 [|T2|. The soliton of ( |1.7[ ) 
has mass 73v/e and size ~ 2/{ev) and we take e as a free parameter since the particles 
corresponding to this soliton have not yet been discovered. 

The standard electroweak Higgs plus gauge boson sector is obtained by gauging the 



SU{2)l X U{1)y subgroup of SU{2)l x SUi2)R in the Lagrangian ([TTl| ). Throughout this 
paper we neglect the U{1) interactions. The complete Lagrangian we consider is obtained 
upon gauging the SU{2)l symmetry of (p!T7|) : 

,2 . _ 1 



1 V 

r = Trpt^" F A Tr 

2c/2 4 



32e2 



Tr 



where 



D^U ={d^-iA,)U 



:i.9) 
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with Afj^ = y4^(T°/2 where the a" are the Pauh matrices. In the unitary gauge, U = 1, and 
the Lagrangian is 

where we have introduced 



m = y and ^ = ^ . (1.11) 



Note that the equations of motion derived from (|1.1CI| ) agree with those obtained by varying 



^ and then setting U = 1. Also note that for fixed m and ^ the classical equations of 
motion are independent of g. Since m is dimensionful and sets the scale, characteristics of 
the classical theory depend only on the single dimensionless parameter i^. 

B. The Soliton and the Sphaleron 

The classical lowest energy configuration of (|1.10|) has A^ = and the quantum theory 
built upon this configuration has three spin-one bosons of equal mass m. In the limit where 
g goes to zero with e and v fixed (hence, ^ goes to infinity) the Lagrangian (|1.8|) is well 
approximated by its ungauged version ( p..7|) which supports a stable soliton, so one suspects 
that for large ^ the Lagrangian ( |1.8|) and its gauge-fixed equivalent ( p..lO|) also support 



a soliton. In fact, Ambjorn and Rubakov showed that for ^ larger than ^* = 10.35 
the Lagrangian ( |1.10|) does support a classically stable soliton whereas for ^ < ^* such a 
configuration is unstable. Let f/i(x) be the winding number one soliton, the skyrmion, 
associated with the ungauged Lagrangian (PTTp. For large ^, this configuration is a good 
approximation to the soliton of the gauged Lagrangian (|1.8|), so in the unitary gauge the 
soliton is A^°^ ~ iUldiUi, Aq°^ = 0. For all ^ > ^* the quantum version of the theory 
described by ( |1.10| ) has, in addition to the three equal mass VT-bosons, a tower of particles 
which arise as quantum excitations about the soliton, just as the proton, neutron and delta 
can be viewed as quantum excitations about the original skyrmion p]T|JT^ . 



5 



The Lagrangian ( |1.10|) determines a potential energy functional which depends on the 



configuration y4^(x). The absolute minimum of the energy functional is at = 0. For 
^ > ^* there is a local minimum at the soliton = A^°^{x.) with nonzero energy given by 
the soliton mass Mgoi- (Of course a translation or rotation of y4^°^(x) produces a configuration 
with the same energy so we imagine identifying these configurations so that the soliton can 
be viewed as a single point in configuration space.) Consider a path in configuration space 
from = to A^°^(x). The energy functional along this path has a maximum which is 
greater than the soliton mass. As we vary the path, the maximum varies, and the minimum 
over all paths of this maximum is a static unstable solution to the classical equations of 
motion which we call the sphaleron of this theory. (The sphaleron of the standard model 
marks the lowest point on the barrier separating vacua with different winding numbers. 
Here, the sphaleron barrier separates the vacuum from a soliton with nonzero energy.) For 
fixed V and e the sphaleron mass Mgph goes to infinity as g goes to zero refiecting the fact that 
for g = 0, configurations of different winding (f/'s with different winding in ( |1.7|) ) cannot be 
continuously deformed into each other. For fixed g and m, as C, approaches ^* from above 
the sphaleron mass comes down until at ^ = ^* the sphaleron and soliton have equal mass. 
For ^ < ^* the local minimum at nonzero energy has disappeared. 



For ^ > the classically stable soliton can decay by barrier penetration [§-[T0|. This 
process has been studied in detail by Rubakov, Stern and Tinyakov who computed the 
action of the Euclidean space solution associated with the tunnelling. They show that in 
the semi-classical limit as ^ — > oo the action approaches 8n'^/g^ whereas as ^ ^ ^* with g 
fixed the action goes to zero since the barrier disappears. 

C. Over the Barrier 

In this paper, we focus on processes where there is enough energy to go over the barrier. 
In the standard model, the sphaleron mass is of order m/g'^ and the sphaleron size is of 
order 1/m. This means that for small g two incident W bosons each with energy half the 
sphaleron mass have wavelengths much shorter than the sphaleron size. This mismatch is the 



reason that over the barrier processes are generally believed to be exponentially suppressed 
in W — W collisions. In contrast, in the model we consider we can take a soliton as one of 
the initial state particles. To the extent that the soliton is close to the sphaleron we have a 
head start in going over the barrier. We can also choose parameters such that an incident W 
boson with enough energy to kick the soliton over the barrier has a wavelength comparable 
to both the soliton and sphaleron sizes. 

We first look at solutions to the Minkowski space classical equations of motion derived 
from the Lagrangian ( |1.10|) . To simplify the calculations we work in the spatial spherical 



ansatz . We solve the equations numerically. As initial data we take a single electroweak 
soliton at rest with a spherical pulse of gauge field, localized at a radius much greater than 
the soliton size, moving inward toward the soliton. In the next section, we display one 
example of a soliton-destroying pulse in detail. For ^ within about a factor of two of for 
all the pulse profiles we have tried with the pulse width comparable to the soliton size, there 
is a threshold pulse energy above which the soliton is destroyed. The energy threshold is 
larger than the barrier height, and does depend on the pulse profile. However, the existence 
of a threshold energy above which the soliton is destroyed seems robust, and in this sense 
the choice of a particular pulse profile is not important. 

A classical wave narrowly peaked at frequency u with total energy E can be viewed 
as containing E/hu particles. Making a mode decomposition, we can then estimate the 
number of gauge field quanta, that is W bosons, in a pulse which destroys the electroweak 
soliton. From (|1.10|) we see that such a pulse has an energy proportional to l/g^ for fixed m 
and ^. Thus, the particle number of any such pulse goes like some constant over g^. For 
example, at ^ = 12 we have found pulses with gf^A^ ~ 2.5. At this value of ^, by varying the 
pulse shape we could reduce g'^N somewhat but we doubt that we could make it arbitrarily 
small. Upon reducing C, towards C,* and thus lowering the energy barrier AE, smaller values 
of g'^N become possible. For example, at ^ = 11 we have found pulses with g'^N ~ 1. In the 
standard model, finding gauge boson pulses which traverse the sphaleron barrier and which 



have small g'^N appears to be much more challenging JTsl. Note from the form of ( |1.10|) 



that taking g to zero with m and ^ fixed is the semi-classical limit. In this limit, the soliton 
mass, the sphaleron mass and their difference AE all grow as l/g"^. The number of particles 
in any classical pulse which destroys the soliton also grows as 1/g'^. 

The existence of soliton destroying classical pulses has quantum implications beyond a 
naive estimate of the number of particles associated with a classical wave. In Appendix B 
we give a full and self-contained account of the relationship between classical solutions and 
the quantum tree approximation in a simple scalar theory. In a theory with an absolutely 
stable soliton, the Hilbert space of the quantized theory separates into sectors with a fixed 
number of solitons, and states in different sectors have zero overlap [|16| . We argue in Section 
III using the results of Appendix B that the existence of classical solutions in which solitons 
are destroyed demonstrates that there are states in the zero and one soliton sectors of the 
quantum theory whose overlap in the semi-classical limit is not exponentially small. These 
states are coherent states with mean number of H^-bosons of order l/g"^. Knowing that some 
quantum processes exist which connect the zero and one soliton sectors suggests that we go 
beyond the semi-classical limit and look for such processes involving only a single incident 
l^-boson. 

There is an interesting limit in which we can reliably estimate amplitudes for single 
particle induced decays. Recall that for m and g fixed, as ^ approaches ^* from above the 
sphaleron mass approaches the soliton mass. We can hold m fixed and pick ^ to be a function 
of g chosen so that as g goes to zero ^ approaches ^* in such a way that AE = Mgph — Mg^i 
remains fixed. We call this the fixed AE limit. It is different from the semi-classical limit 
in that as g goes to zero the classical theory is changing. We will argue in Section IV that 
for ^ near ^* it is possible to isolate a mode of oscillation about the soliton whose frequency 
is near zero, which is in the direction of the sphaleron. This normalizable mode, which 
we call the A-mode, is coupled to a continuum of modes with frequencies u > m. If the 
A-mode is sufficiently excited by energy transferred from the continuum modes, then the 
soliton will decay. We are able to estimate the amplitude for a single VT-boson of energy 
E to excite the A-mode enough to induce the decay. At threshold the cross section goes 
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like exp(— c/(7^/^) where c is a dimensionless constant. In the same hmit we can calculate 
the rate for the soliton to decay by tunnelling and we get exp(— (9/(9 — 2-\/3))c/(7^/^). Both 
are exponentially small as g goes to zero and the ratio of the tunnelling rate to the induced 
decay rate is exponentially small. 

D. Fermion Production 

We introduce fermions into this theory as in the standard model. The left-handed com- 
ponents transform as SU{2)l doublets whereas the right-handed components are singlets. 
The fermion mass is generated in a gauge-invariant way by a Yukawa coupling to the Higgs 
field. For simplicity we only consider the case where both the up and down components of 
the fermion doublet have equal mass mj. In any process where a soliton is destroyed there is 
a violation of fermion number. The nature of this violation is different depending on whether 
the fermion is light, mfL ^ 1, or heavy, mfL ^ 1, where L is the characteristic size of 
the soliton. In the light fermion case, when the soliton disappears one net anti-fermion is 
produced in the process. In the heavy fermion case no fermions are produced. However in 
this case the soliton carries heavy fermion number and when the soliton is destroyed this 
quantum number is violated. In both cases there is a change of fermion number of minus 
one and heavy minus light fermion number is conserved as it must be since the heavy and 
light fermion number currents have the same anomalous divergence. 

E. Relating the Model to the Real World 

The metastable electroweak soliton of the modified Higgs sector is an intriguing object 
to study. Yet this beast is not found in the standard electroweak theory where the Higgs 
sector is a linear sigma model with no higher derivative terms. It is reasonable to ask if 
the modified theory gives a credible description of physics at the weak scale. To date the 
Higgs boson has not been found. If it is found and the mass is low so that A of (|1 . 1|) is 
small then working in the infinite A limit would not well approximate reality. However, if 
the Higgs is heavy, then working with infinite A could be justifiable. Working at the scale 
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V and below, we then integrate out the heavy Higgs, leaving a low energy effective action. 
In this strongly interacting case, higher derivative terms in the effective action would not 
be perturbatively small and we would expect all possible higher derivative terms consistent 
with the symmetries. This effective theory would or would not support stable solitons. If it 
did then our use of the Skyrme term is justified as a simple way to write an effective action 
which supports solitons. 

It is possible that the Higgs is not fundamental. Rather the Higgs sector may be an 
effective theory describing the massless degrees of freedom which arise as a result of spon- 
taneous symmetry breaking in some more fundamental theory. For example, this is the 
basis of technicolor theories in which the symmetry breaking is introduced via a scaled up 
version of QCD. In technicolor theories one finds technibaryons which can be described as 
electroweak solitons just as the baryons of QCD can be described as skyrmions. For now, 
regardless of whether the underlying theory is specifically a technicolor model, as long as 
we are consistent with symmetry considerations, we are free to choose the effective theory 
to conveniently describe the particles which interest us. Thus ( |1.7| ) is a simple way to de- 
scribe three massless bosons (which are eaten in the gauged version ( |1.8[ )) as well as a stable 
(metastable in (|1.8| )) heavy particle. Of course the effective theory includes higher derivative 
terms other than the Skyrme term, so it is not the precise form of ( |1.8|) which we think is 
plausible, but rather the physical picture which it describes. 

It is worth asking what processes can sensibly be described using the effective theory. 
The effective theory is a derivative expansion in momenta over v. Consider the (fermion 
number conserving) production of soliton - anti-soliton pairs in ly — ly collisions. These 
processes are beyond the regime of applicability of the effective theory because the incident 
particles have momenta which are greater than f , and the underlying theory must therefore 
be invoked. (For example, in a technicolor theory the production process would be described 
as techni quark - anti-techni quark pair production followed by technihadronization.) The 
effective theory is, however, well-suited to describing soliton decay induced by a single W 
boson with energy just above AE in the fixed /S.E limit. In this limit m is held fixed while 
5^ — >■ 0, and thus v —>■ oo. Therefore, the ratio of the incident W momentum to the scale v 
is going to zero, and a treatment using the effective theory is justified. 
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II. SOLITON DESTRUCTION SEEN IN CLASSICAL SOLUTIONS 



We begin our investigations classically. We wish to find solutions to the Minkowski 
space classical equations of motion derived from the Lagrangian ( |1.10|) which at early times 
have an electroweak soliton and an incident pulse and which at late times have outgoing 
waves only, the soliton having been destroyed. In this section, we investigate solutions to 
the equations of motion numerically. In order to make the numerical problem tractable, we 
work in the spatial spherical ansatz [p^ . 

The unitary gauge Lagrangian ( |1.10| ) yields the equations of motion 

1 



(2.1) 



where 



D^F>''' = d,F>'''-t[A^,F>'''] . 



(2.2) 



In the unitary gauge Lagrangian, the Skyrme term, Tr [A^, A,^] , is the same as the quartic 

^2 



term in TrF^^^. Thus the unitary gauge equations of motion ( |2.1| ) are the same as the 



equations of motion for a massive non-abelian vector field except that the coefficient of the 
cubic term is now (1 + 1/4,^). The classical equations of motion depend only on m, which 
sets the scale, and on the dimensionless parameter ^, but do not depend on g. 



The spherical ansatz [n] is given by expressing the gauge field in terms of four real 
functions Oq , ai , a and 7 of r and t: 



A(x,t) = -ao(r,t)a- X 



A(x,t) = - 



ai(r, tjcr ■ xxj H (cij 



^(r, t) 

a ■ 9.Xi) H — eijkXjak 



(2.3) 



where x is the unit three-vector in the radial direction and a are the Pauli matrices. For A^ 
to be regular at the origin, we require that oq, a, ai — a/r, and 7/r vanish as r ^ 0. In the 

are 



spherical ansatz, the unitary gauge equations of motion ([O 



00 -L . 

r m H Y 

2r 



[2.4a) 
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i\x\ 



X = -m^ (X + + ^ (X + 



(2.4b) 



where 



X = « + n7 



(2.5a) 
(2.5b) 



Df,x = {df,-ia^)X ■ (2.5c) 
The indices take the values and 1 and are raised and lowered with the 1 + 1 dimensional 

metric ds"^ = dt^ — dr"^. The notation suggests that we are dealing with a 1 + 1 dimensional 

U{1) gauge theory with gauge field and a complex scalar x of charge 1. In fact the left 

hand sides of ( p.4|) are U{1) gauge covariant whereas the right hand sides involving the mass 

and Skyrme terms are not. This can be understood as follows. Before gauge fixing, the 

underlying theory ( |1.8| ) is SU{2) gauge invariant. If we take fields in the spherical ansatz, 

gauge transformations of the form exp [iQ{r, t)a ■ x/2] keep the fields in the spherical ansatz. 

Thus, the spherical ansatz has a residual U{1) gauge invariance. In the unitary gauge the 

mass and Skyrme terms lose their covariant form which is seen in ( p. 41) . Note that ao, which 

determines Ao, is not a dynamical degree of freedom so the problem has been reduced to 

the dynamical degrees of freedom x and ai. 

Our task is to choose initial conditions and then to evolve the fields forward in time. We 
specify initial conditions by specifying x and ai and their time derivatives at t = 0. We 
then use the z/ = component of ( f^.4aD which is Gauss' law to determine Oq. The u = 1 
component of ( |2.4a| ) and the equation ( |2.4bD are the second order equations of motion which 
we use to evolve x aiid ai forward in time. At every time step, we update oq using Gauss' 
law. We present some details of the numerical methods in Appendix A. In order to have a 
check on the accuracy of our numerical methods we have also solved the equations in Aq = 
gauge and the details of this approach are also given in Appendix A. 

In describing the solutions, it is convenient to write x as 



X = -ip{r,t)exp [i9{r,t)] 



(2.6) 
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The r = boundary conditions on a and 7 imply 

p(0,t) = l (2.7a) 

^(0,t) = 0. (2.7b) 

The boundary condition ( p. 7b ) should strictly be that 6'(0,t) is an integer multiple of 27r. 



However, since p never vanishes at the origin, Q is constant in time there and we have taken 
it to vanish. In vacuum, p = 1 and ^ = everywhere. Finite energy solutions must satisfy 

lim 0(r,t) = 2n7r (2.8) 

at all times. Thus we see that in the spherical ansatz finite energy configurations with p 7^ 
at all r can be characterized by n, the integer-valued winding of the x field. This winding is 
the number of times the complex-valued x wraps around x = as x goes from — i at r = 
to — i at r = 00. Note that this winding can change only if x passes through x = 0? that is, 
if p goes through zero at some t and r. 

We now look at the soliton in terms of the variables x aiid a^. To understand the 
qualitative form of the soliton configuration, it is useful to begin with the Skyrme model as 
we did in Section I. Recall that for (7 ^ with e and v fixed, the Lagrangian (|1.1CI| ) reduces 
to ( |1.7| ) and the soliton becomes the Skyrme soliton written in unitary gauge. In this limit, 
^soi _^ i\j\di\J\ (v4q°^ = for all values of g) where V\ is the winding number one Skyrme 
configuration. Now \J\ is of the form 

f/i =exp[io'-xF(r)/2] (2.9) 

where -F(O) = and F{oo) = 2tt. In terms of x and this configuration is 

a/. = , x = -^exp[^F(r)] . (2.10) 

In this case we see that n, the winding of x, is equivalent to the winding of f/, both of which 
are 1. For (7^0 with e and v fixed, the electroweak soliton configuration is described by 
( p.lOp . For nonzero g with > the soliton field configuration is still approximately of 
the form (|2l0| ). 
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To find the precise form of the electroweak soliton we add energy non- conserving damping 
terms to the equations of motion. Specifically, we add —Tdx/dt and —Tdai/dt to the right- 
hand sides of ( p.4b| ) and the u = 1 component of ( |2.4a| ) respectively where F is a constant. 



The solutions to the modified equations of motion lose energy as they evolve. Depending 
on the initial configuration, this modified evolution leads either to the vacuum or to the 
soliton.Q For a given ^, we find the soliton by choosing initial configurations with n = 1 and 
evolving them using the modified equations. When we find an initial configuration which 
evolves to a nonvacuum configuration, we check that the configuration so obtained is indeed 
a static stable solution to the unmodified equations of motion. In Figure 1, we show p, 6, 
X and the energy density for the electroweak soliton with ^ = 12. In Figure 2, we show 
how Xsoi changes with ^. For ^ approaching the soliton configuration does not change 
qualitatively, and in particular p remains well away from zero and n remains 1. In the 
^ — s> oo limit, p ^ 1, ai — s> like m/y/^, and the size of the soliton, i.e. the size of the 
region over which 6 varies, shrinks like l/{my^). 

We now consider initial conditions with a soliton and an incoming pulse which destroys 
the soliton. We have experimented with several ansatze for the pulse shape. Here we present 
one which we feel is fairly simple and which does destroy the soliton. Recall that for the 
soliton X = — z at r = and wraps once around the origin as r increases from to infinity. 
At t = 0, the incident pulse we choose has Xpuise = Xvac = —i for r < (where tq is large 
compared to the soliton) and has Xpuise Xvac as r ^ oo. As r increases from r^, we choose 
Xpuise to loop in the complex plane in the opposite direction to that in which Xsoi winds. 
Xpuisc is a small enough excitation about Xvac that Ixpuise — XvacI < 1 and the pulse has n = 0. 
Specifically, as the initial condition for x we use 

X = (Xsol - Xvac) + Xpulse (2.11) 



^It can also lead to a multiple winding number soliton with n > 1. These configurations have 
been studied by Brihaye and Kunz They have more energy than n widely separated solitons, 
and we are not concerned with them in this paper. 
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FIG. 1. The soliton configuration for ^ = 12. Panels a) and b) show p and 6 as functions 
of r measured in units of the inverse VF-mass. Note that changes by 27r. Panel c) shows 
X = —ipexp{i9) vs. r. The projection of this curve onto the x-plane is shown in Figure 2. In 
panel d) we show the energy density vs. r. By energy density we mean times the 3-dimensional 
energy density; the total energy of the configuration is the area under this curve. 
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-1 -0.5 0.5 1 

Rex(r) 

FIG. 2. X foi' the electroweak soliton for different values of ^. All the curves begin at x = ~* 
at r = 0, traverse a counter-clockwise loop which encircles x = 0, and return to x = ~^ ^-s goes 
to infinity. The dotted, dot-dashed, dashed, and solid lines correspond to = oo, 15, 12, and 10.5 
respectively. Recall that there is no stable soliton for ^ < ^* = 10.35. The range of r over which 
the loops are traversed, i.e. the size of the soliton, is approximately 2/(m^^/^). At ^ = oo, |x| = 1 
and the dotted curve is a circle. This reflects the fact that as ^ approaches infinity the soliton 
configuration approaches the form (2.10). 

where 

Xpuise = + I [e^-^^^'-) - l] , (2.12) 
with h a constant whose absolute value is less than one and with the function g{r) given by 

(^) = |exp[-(r-ro)2M ; r > ^2.13) 
I 1 ; r <rQ 

We choose initial conditions with ai = a^°' and with ai such that = 0. Now, we must 
specify x- We wish to do this in such a way that the energy of the pulse propagates inward 
toward the soliton rather than outward toward large r. In a massive theory, it is impossible 
to achieve this exactly, and we use the following prescription which works well enough for 
our purposes. For x = Xpuise given by ( p.l2[ ) and ( |2.13| ), (x — Xvac) has a mean wave number 



squared of approximately (tt/ct)^ and a mean frequency u ~ ym^ -|- {n/a-y. So, we define 
a velocity 
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and choose the initial condition 

X = v^^. (2.15) 

Thus, in this ansatz the initial conditions are parametrized by the amplitude b, the pulse 
width (7 and the initial radius tq. 

Figure 3 shows the result of hitting the sohton at ^ = 12 with a pulse chosen using the 
above ansatz with b = 0.23, a = 7r/5m, and tq = 4/m. We plot p, 6, x and the energy 
density as functions of r for a number of different times. First, note that the pulse does 
move inward toward the soliton. The soliton energy is 72.67 m, the sphaleron energy is 
73.9 m, and the total energy of the solution is 85.85 m. Thus, neglecting the small amount 
of energy initially in the pulse which goes outward, the soliton is hit with a pulse with 
energy -Bpuise = 13.18 m which is larger than AE — Mgph — Mgoi — 1.2 m. We see that at 
time t = 7.2/m there is an outgoing pulse with somewhat less energy than -Epuiscj and the 
soliton has been somewhat distorted, as it begins to fall apart. At time t — 8.59/m, p is 
very close to zero at r = 0.94/m. At late times, we see from the plot of 9 that the winding 
n is zero, and we see from the plot of the energy that there is in fact no soliton present. We 
estimate N, the number of FT-bosons in the incident pulse which destroyed the soliton, as 

g^N ~ ~ , ^P"'^*^ ~ 2.5 . (2.16) 

We now sketch how the results of Figure 3 change as we vary b and a. First, upon 
increasing b from 0.23 (and thus increasing E'puise and g^N) the time delay between the 
emergence of an outgoing pulse and the collapse of the soliton decreases - the soliton is 
destroyed more promptly. Upon decreasing 6, the time delay increases — less energy is 
dehvered to the soliton and the soliton takes longer to fall apart. Decreasing b further, we 
find a threshold somewhere between b = 0.21 and b = 0.22 below which the soliton survives. 
Below this threshold, after the emergence of an outgoing pulse, the soliton radiates any 
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FIG. 3. This figure, which continues on the following pages, shows the destruction of the ^ = 12 
soliton by the pulse specified in the text. The simulations were done on a lattice which extended 
to r = 16/m, but only < r < 12/m is displayed. On this page, we show snapshots of p{r) at 
eight different times. Each panel displays p{r) at two times; the solid curve shows p{r) at the 
earlier time, and the dot-dashed curve shows p{r) at the later time. For example, in the first panel 
the solid curve at t = 0.0/m shows the incident pulse superimposed on the soliton of Figure 1. 
The dot-dashed curve shows p at t = 3.2/m when the pulse has moved inward toward r = 0. By 
t = 4.8/m, the solid curve in the second panel, it is clear that the soliton has been disturbed. At 
t = 7.2/m, there is an outgoing pulse at r ~ 3/m and the soliton has not returned to its initial 
shape. At t = 8.59/m, p = at r = 0.94/m. At t = 15.2/m, the final time shown, there is an 
outgoing pulse at r ~ 11/m followed by the outgoing remnants of the soliton at r 6/m. 
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t = 0.0, 3.2rTV,"^ t = 4.8, 6.4rTVv-i 




FIG. 3. On this page, wc show snapshots of 9{r) corresponding to the preceding snapshots of 
p{r). We show two times per panel, and the sohd curve is the eaxher of the two. At all times before 
t = 8.59/m, 9 increases from ^ = at r = to ^ = 27r at large r. At t = 8.59/m, note the large 
slope in 9 dX r = 0.94/m which is where p vanishes. At later times, 9 is no longer wound. 
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t = 0.0 



t = 3.2 mw"^ 




FIG. 3. On this page and the next, we combine p and 9 from the previous pages into 
three-dimensional plots of x = —ipexp{i9) as a function of r at eight times. Initially, we see 
the pulse incident upon the soliton. Between t = 6.4/m and t = 9.6/m, we see the soliton shed 
an outgoing pulse, and then shrink from a loop which encircles x = to an excitation about the 
vacuum x = ~^ which does not. Att = 15.2/m,xis close to the vacuum for small r, and at larger 
r we see the outgoing pulse and the outgoing remnants of the soliton. 
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21 



t = 0.0, 3.2 171^-^ 




4 6 8 

radius (rriw ) 



12 



200 



150 



100 



t = 4.8, 6.4m,;v" 




2 4 6 8 10 12 



t = 7.2, 8.59 m^v- 




200 r 



150 



100 



t = 9.6, 15.2rT\v" 



2 4 6 8 10 12 




/ \ 



8 10 12 



FIG. 3. Finally, we show snapshots of the energy density at the same times as before. There 
are two times per panel, the solid line at the earlier time and the dot-dashed line at the later time. 
As in Figure 1, by energy density we mean times the 3-dimensional energy density. At t = 0, 
we see that the incident pulse has much less energy than the soliton, but nevertheless it destroys 
the soliton. The destruction of the soliton is seen most clearly by looking at t = 15.2/m. At this 
time, there is no energy density visible for r less than about 3/m, which is where the soliton was 
at t = 0. The soliton is no more. There is an outgoing pulse at r ^ 11 /m which has an energy 
comparable to but slightly lower than that of the incident pulse. This is followed at r ~ 6/m by 
an outgoing shell with energy comparable to that of the soliton. 
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remaining excess energy outward and settles back to its undisturbed state. For several 
values of a ranging from half to twice that in Figure 3, the threshold h is about the same. 

We have worked at values of ^ ranging from 10.5 to 100. For a given a the threshold 
amplitude is lower for values of ^ closer to For ,^ = 11, for example, we have found soliton 
destroying pulses with g'^N ~ 1. As ^ becomes very large the soliton size (~ l/{my/^)) 
becomes much smaller than the sphaleron size (~ 1/m) and the barrier height AE grows 
like Msoiv^. At large ^, therefore, the energy of soliton destroying pulses must become large 
compared to Mgoi and large compared to the inverse sphaleron size. It is nevertheless a 
logical possibility that such pulses could be found with high frequencies and small values of 
g'^N . For ^ = 50 and above, however, we have only found soliton destroying pulses which 
have large g'^N . This suggests that because at large ^ the soliton is no longer similar to the 
sphaleron, we lose the advantage that we have in this model, relative to the standard model, 
in finding sphaleron crossing solutions with g'^N small. 

We have chosen to present results for ,^ = 12 (rather than choosing closer to ^* where the 
threshold amplitude and threshold g'^N are lower) because at ^ = 12 the tunnelling lifetime 
of the soliton is much longer than any time scale in Figure 3. As we discuss in Section V, 
Rubakov, Stern and Tinyakov |T^ write the tunnelling lifetime as as r ~ (l/m)exp(2i?) 
and calculate that g^B = 4 ± 1 for ^ = 12. Thus, for g = 0.65 and ^ = 12, r ~ 10^/m. 

We have tried a number of incident pulse shapes that do not fall into the ansatz we have 
described in detail, and have found qualitatively similar results. For all the cases which we 
have considered with ^ within a factor of two of we have observed that as we vary the 
amplitude of a pulse whose size is comparable to that of the soliton, for amplitudes above 
some threshold the soliton is destroyed. For a given pulse shape, the threshold energy and 
threshold g'^N decrease as ^ decreases toward C,*- Among the few pulse shapes which we 
tried, the threshold energy was lowest for the ansatz of ( p.l2| ), but we have certainly not 
found the lowest energy or lowest particle number pulses which destroy the soliton. Indeed, 
a soliton destroying pulse with energy just above AE could be obtained by starting with a 
slightly perturbed sphaleron, watching it decay to the soliton, and then time reversing. We 
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will see in Section IV that for ^ near ^* a "pulse" so obtained would be a very long train 
of small amplitude waves, rather than a simple pulse of the kind we have used to destroy 
solitons in our numerical experiments. 

The lesson of this section is that in the model we are treating, it is straightforward 
to find soliton destroying, sphaleron crossing, fermion number violating classical solutions. 
Particular pulse profiles are not required — pulses of any shape we have tried (with sizes 
comparable to the soliton size) destroy the soliton if their energy is above some shape- 
dependent threshold. 

III. QUANTUM IMPLICATIONS OF CLASSICAL SOLUTIONS WHICH 

DESTROY SOLITONS 

In a theory like the ungauged Skyrme model with a static classical solution which is 
absolutely stable, that is, separated by an infinite energy barrier from the classical vacuum 
configuration, the Hilbert space of the quantized theory separates into sectors with a fixed 
number of solitons, and states in different sectors have zero overlap [|T6|. The one soliton 
sector, for example, is a Fock space of states with one soliton and any number of mesons. 
The mesons (pions in the Skyrme model) are the quantized fluctuations about the soliton 
configuration, and the states in the one soliton sector are scattering states of mesons in the 
presence of a soliton. No process, not even one involving large numbers of mesons, connects 
states in the one soliton sector with states in the vacuum sector. 

In our theory, the electroweak soliton is not absolutely stable. It is separated by a 
barrier of finite height from the vacuum. The Hilbert space has sectors with a fixed number 
of solitons and any number of VT-bosons. However, we now argue that the existence of 
the classical solutions described in the previous section, in which incident pulses destroy a 
soliton, demonstrates that there are states in the zero and one soliton sectors with nonzero 
overlap. 

Consider a classical solution obtained by taking a solution in which a soliton is destroyed 
and the time reverse of a different such solution and combining them as we now describe. 
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At very early times, there are two incoming pulses widely separated in time. The inner 
pulse, of total energy Ei, is the time reverse of a solution of the kind found in Section II. It 
forms a soliton, an outgoing pulse of energy Ei — Mgoi is radiated, and the soliton of mass 
Msoi is left sitting at the origin. Then the outgoing pulse passes the second incoming pulse 
at a radius large enough that the amplitudes of both pulses are small, and no interaction 
occurs. Subsequently, the second pulse of total energy E2 arrives at the soliton and destroys 
it, yielding a second outgoing pulse of energy E2 + Msoi. At very late times there is no 
soliton present, and there are two outgoing pulses. This entire solution falls into the class 
of classical solutions discussed in Appendix B, in that at very early and very late times the 
fields are small amplitude excitations about the vacuum. By the arguments of Appendix B, 
the existence of this solution implies that we can construct normalized coherent states 

\fLf-n) (3.1) 

and 

I/out' /out) (3-2) 

such that 

(/out, /outl/in, fl) = + ^(^°)} (3-3) 

as (7 ^ 0, with 6 a real phase. The energy and particle number in the in and out states are 
0(1/(7^). Equation (|B24| ) expresses 6 as an integral of the classical fields over space-time. 
When the time separation T between the pulses 1 and 2 is large compared to the times 
necessary for the formation and destruction of the soliton, this integral can be written as 
the sum of three terms, 

^ = 4 + 4 + ^soiT 3.4 

r 9^ r 

where 61 depends only on fl^^^ and the formation solution, 62 only on f?^ and the destruction 
solution, and Msoi is the classical energy of the soliton which is of order l/g"^. 
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With this information concerning asymptotic states, the only possible interpretation is 
that there are coherent states of H^-bosons in the one soliton sector |sol, /^^J and |sol, /j^), 
and that 



(sol,/iJ/i)=exp|4 + 0(^' 



(3.5) 



9' 



(3.6) 



Thus there are processes connecting the one-soliton sector to the vacuum sector which are 
not exponentially suppressed as g ^ 0, though they do involve O^l/g"^) W^-bosons. 

In the remainder of this paper our goal is to study quantum processes in which a single 
W^-boson incident upon the soliton kicks it over the barrier causing it to decay. In Section 
V, we will do a controlled calculation of this process in a limit in which ^ goes to C,* as g 
goes to zero. In order to do this calculation, however, we first need a better understanding 
of classical dynamics of the theory with ^ near and to this we now turn. 

IV. CLASSICAL DYNAMICS FOR ^ NEAR 

In order to discuss the special features of the dynamics of our system for ^ near and 
because we will need it to discuss the quantum version of this theory, we introduce the 



Hamiltonian which arises from (|1.10|) 



H 



d^x <^ g^Ti WW + — 
I 9^ 



Ixr F'^F'' - m'Ti A^A^ - ^Tr [A^ , A,f 



2Tr 



(4.1) 



where 



An* = d.W -i Ai,W 



(4.2) 



Now A^ has no conjugate momentum and the A^ equation is 
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+ ^[[A\ A^] , A,] + g^DiW = . (4.3) 

This linear equation for A° can be solved giving y4° in terms of A* and 11* but we do not need 
to do this explictly. The Hamiltonian for our system is given by ( [4.1|) with A^ determined 
by ( [4. 3D and has the general form 

H = ^UM-\A)U+^V[A] , (4.4) 

where the sum over the coordinate x, the spatial index i and the group index are all implicit. 
The matrix M~^{A) involves derivatives with respect to x, depends on the configuration A, 
and we assume that M'^i^A) is positive. Note that static solutions to the equations of 
motion, that is those with U = A = 0, occur where SV/6A = and have 11 = 0. The 
classical equation of motion for A which arises from (|4.4| ) is independent of g. Thus for the 
discussion of classical dynamics which we are having in this chapter, we can set g = 1. We 
will restore the g dependence in the next chapter. 

The potential energy functional V[A] has its overall scale set by m but the topography of 
fixed energy contours is set by ^. Ambjorn and Rubakov |]10| showed that for ^ > = 10.35 



there is a local minimum, the soliton, whereas for ^ < ^* this minimum is absent. For ^ > ^* 
there is also a sphaleron, that is a saddle point configuration whose energy is greater than 
that of the soliton. As ^ approaches ^* from above, the sphaleron and soliton merge. 

We are particularly interested in configurations which, at least initially, are small pertur- 
bations around the soliton. To work with these configurations we find it convenient to make 
a canonical transformation which has the effect of setting M~^(y4soi) = 1 and 



= 0. 

^sol 



dA 

To see that this is possible let fa be some complete set of orthonormal, spatial vector, 
matrix- valued functions of x, indexed by a, which can be used to expand 11 and A. Let 
the coefficients of the expansion of A relative to the soliton be and the coefficients of the 
expansion of 11 be that is 

A(x,t)-A-^(x) = ^g"(t)/,(x) 

a 

n{^,t) = J2pait)u^) . (4.5) 

a 
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(Note that the transformation from A{'K,t), n(x, t) to q"'{t), Pa{t) is canonicaL) Upon 
making this transformation, ([4 .41) has the form 



H=^g''^iq)PaPp + Viq) . (4.6) 
A canonical transformation of the form 

r = q"'{q) and P'a = ^Pp (4-7) 

can be viewed as a general coordinate transformation with pa transforming as a covariant 
vector. It is always possible to choose coordinates such that 

is equal to 5"^ with dg"^^ /dq'^ = at any given point. In fact this can be accomplished at 
= (the soliton) with a transformation of the form q'" = C'^q^ + Cp^q^q^ . This means 
that the Hamiltonian ( |4.6| ) can be written as 

H =]^p^[5^P + 0{q')\pp + V{q) , (4.9) 

where we have made the required canonical transformation and dropped the primes. Note 
that V{q = 0) = Msoi and ((9r/9g°)|g=o = 0. 

For ^ > ^* consider small oscillations about the soliton. The frequencies squared are given 
by the eigenvalues of the fluctuation matrix d'^V / dq'^dq^ at g = 0. The soliton is a localized 
object so fluctuations far from the soliton propagate freely. Therefore the fluctuation matrix 
at the soliton has a continuous spectrum above m?. A given soliton configuration and a 
translation or rotation of that configuration have the same energy and both solve dV/ dq"' = 
0. This implies that at g = there are six zero eigenvalues of d'^V/dq°'dq^ . The associated 
modes which correspond to translating and rotating the soliton are not of interest to us and 
will be systematically ignored. 

For ^ close to ^* we now argue that there is one normalizable mode whose frequency ujq 
goes to zero as ^ goes to To see this we write 
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dV 



dV 



Qsph 



+ 



q=0 



dq'^dqf^ 



Vsph ^ 



(4.10) 



q=0 



2 dq'^dqPdq' 

At the soliton (g = 0) and at the sphaleron the first derivatives are zero. As ^ approaches ^* 
the sphaleron and sohton merge so g^j^ goes to zero. It is useful to introduce the normalized 
function q^^h 

5'sph 



9'sph 



Q 



(4.11a) 



where 



^sph^sph ■ 



(4.11b) 



As goes to Q goes to zero but gsph does not. From ( |4.1CI| ) we then have 



dq^'dql^ 



-a -P 
'?sph5'sph 



1 93 V 

1 



g=0 



(4.12) 



g=0 



2^ dq'^dq'^dq' 

For ^ > the fluctuation matrix d'^V/dq°'dq^ at the soliton has only positive eigenvalues 
(except for the translation and rotation zero modes which play no role in this discussion). 
Equation (|4.12| ) tells us that at ^ = ^* where Q = 0, the fluctuation matrix has a zero 
eigenvalue with eigenvector ggph whereas for ^ close to ^* there is a small eigenvalue, Uq, 
whose associated eigenvector is close to gsph- Note that gsph points from the soliton to the 
sphaleron. Thus the low frequency mode, which we call the A-mode, is an oscillation about 
the soliton close to the direction of the sphaleron. 

For ,^ > ,^*, at the sphaleron there is one negative mode, that is one negative eigenvalue 
of the appropriately defined fluctuation matrix. As comes down to C,* the sphaleron and 
soliton become the same configuration so this negative eigenvalue must come up to zero 
in order for the spectra of the fluctuation matrices of the soliton and sphaleron to agree 
at ^ = Therefore for ^ close to ^* the unstable direction off the sphaleron has a small 
negative curvature. There are two directions down from the sphaleron. One heads toward 
the soliton and the other heads (ultimately) to the classical vacuum at A = 0. We see that 
for ^ near ^* the soliton can be destroyed by imparting enough energy to the A-mode since 
it is this mode which is pointed towards the sphaleron and beyond. 
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We wish to describe the interaction of the A-mode with the other degrees of freedom. 
We use the Hamiltonian written in the form ( [4.9| ). At this point it is convenient to make an 
orthogonal transformation on the {q""} so that the transformed set are the eigenvectors of 
the sohton fluctuation matrix d'^V/dq°'dq^\q=o- We will label these vectors as q^^ where 
is the eigenvalue of the fluctuation matrix. The eigenfunctions include: 

i) The continuum states q^^ with eigenvalues cj^ > m^. (Note that for each a;^, in 
general, there is more than one eigenvector. The extra labels on q^ are suppressed 
in our compact notation.) 

ii) The normalizable state = A with eigenvalue which goes to zero as ^ goes 
to C- 

iii) The zero eigenvalue states associated with translation and rotation. 

iv) Other normalizable states which might exist but whose frequencies do not have 
any reason to approach zero as ^ goes to 

Up to cubic order the Hamiltonian ( |4.y| ) is 

+ ]- j dupl + l- f dujLj'^ql+ f dudu'du" c{uj,uj',uj")q^q^'q^» 

+ A^ / dud{uj)q^ + \ dudu' e{uj,uj')q^q^f + . . . (4.13) 

Jm Jm 

where in the ellipses we now include all terms with modes of type iii) and iv) as well as higher 
order interactions of the A-mode and the continuum modes, p is the momentum conjugate 
to A and p^) is the momentum conjugate to q^i- The number h and the functions c, d and 
e are determined by the soliton configuration. For example d{(jj) is presumably peaked at 
values of uo which correspond to wavelengths of order the size of the soliton. As ^ goes to 
^* we know that uiq goes to zero but we expect no dramatic behavior of 6, c, d or e in this 
limit. 

Consider the A-mode potential 
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ViX) = lujlX' + lx' + ... . (4.14) 

There is a local minimum at A = which is the soliton and a local maximum at A = — c<Jq/6, 
which is approximately the sphaleron, where the second derivative is —uJq. We work with ^ 
sufficiently close to ^* that ujq is small. This means that A at the sphaleron is small and if we 
only study dynamics up to and just beyond the sphaleron we are justified in neglecting the 
quartic and higher terms in A. We also see that as ^ goes to ^* so that Uq goes to zero, the 
soliton and sphaleron come together and at ^ = ^* the A potential has an infiection point at 
A = and the soliton is no longer classically stable. 

In order to discover the relationship between uq and — C,*) as ^ approaches it is 
necessary to study the behavior of the A- mode potential as ^ approaches In ( [1.14| ) for 
every value of we have shifted A so that the minimum of the potential is at A = 0. This 
^ dependent change of variables obscures the behavior of the coefficients of the potential 
before the shift. Calling the unshifted variable A, then if we expand the potential in terms 
of e = — ^* about e = where there is an infiection point, we have 

l^(A, e) = 0{e)X + 0{e)X^ +{b + 0{e))X^ + . . . , (4.15) 

where 6 is a constant. We know that the coefficients of A and A^ are zero at e = 0, and we 
assume that these coefficients can be expanded about e = and we know of no reason for 
the order e terms to vanish. For e > the minimum of the potential is at A ~ e^^^, (A is 
shifted relative to A by this amount), and at the minimum of the potential d'^V/dX^ ~ e^/^, 
that is 

^o~(e-r)'/'- (4.16) 

A small amplitude oscillation of the A mode will decay because of its coupling to the 
continuum modes which can carry energy away from the soliton. However for uq < m this 
decay is very slow in the sense that the characteristic time for the decay is much greater than 
I/ujq. To understand this consider A(t) as a source for radiation in the continuum via the 
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coupling dLjd{u!)qi^ in the Hamiltonian ( |4.13| ). Suppose that X(t) is a purely sinusoidal 
oscillation with frequency c^o and with an amplitude which is small. Radiation with frequency 
ujo is not possible because the continuum frequencies begin at = m. However, has 
frequency 2ujo and therefore if cuq > itl/2 the coupling will excite propagating modes with 
uj = 2ujq and the A oscillation will radiate at twice its fundamental frequency. Because the 
coupling is of order A^, the rate of energy loss will be small. \i ujq < m/2 then radiation at 
uo = 2ujq is also not possible. However, if m/3 < Uq < m/2 the A'^g^ coupling (which we have 
not written in ( |4.13|) because it is fourth order) allows the A oscillation to radiate at three 
times its fundamental frequency. There is another source of radiation with uj = Suo. The 
potential for the A-mode is not exactly quadratic so the A oscillation, although periodic, is 
not exactly sinusoidal. If the period of the oscillation is 2tt/ujo, A will be a sum of terms 
of the form sincjot, sin2c(Jot, sin 3co'o''^r-- with diminishing coefficients. This means that A^ 
will also be a sum of terms of this form. Those terms in A^ with frequencies greater than 
m will excite radiation via the X^q^^ coupling. As ujq is reduced from m toward zero, the 
radiation is produced only by higher order couplings and by higher harmonics, and therefore 
the amplitude is reduced and the decay takes longer. 

We have numerical evidence for this behavior within the spherical ansatz. To watch 
an oscillating soliton radiate for a long time, we implement energy absorbing boundary 
conditions at the large r boundary of the simulation lattice, as described in Appendix A. 
We wish to excite the A-mode and watch it oscillate. It is convenient to choose initial 
conditions by starting with some n = 1 configuration and evolving it using the equations 
of motion with damping terms added as described in Section II. Instead of running for long 
enough that the configuration is damped down to the soliton, we stop somewhat earlier. This 
yields a configuration which is the soliton plus a small perturbation. Because the damping 
terms damp modes with higher frequencies more quickly than those with lower frequencies, 
the perturbation that remains is mostly in the lowest few modes. We use the configuration 
just described as the initial condition for the equations of motion with no damping terms. 
The resulting evolution is shown in Figures 4 and 5 for ^ = 10.4. The functions p, 6, and 
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FIG. 4. As described in the text, we have perturbed the ^ = 10.4 sohton and let it evolve for 
a long time. Here, we show p{r) for a series of different times: t = 0, 144, 288, . . . 1440 m~^. This 
shows the envelope of the oscillation of p. In Figure 5, we show p aX r = 0.608/m and r = 10/m 
as a function of time. 

ai (we show p only) all oscillate about the values they take at the soliton and the period 
of oscillation is 16.69 m"^. We identify this with the A-mode and so obtain ujq = 0.3764m. 
Furthermore we see that away from the soliton there is a small amplitude train of outgoing 
radiation. After a brief initial period during which any perturbations not in the A-mode 
radiate away, the outgoing radiation settles down to a frequency 1.129 m, three times the 
fundamental frequency. (At r = 10 m~^, we see in Figure 5 that the frequency 3u;o oscillation 
of p has a small modulation with frequency ujq. This is the tail of the A-mode oscillation 
and is not seen at larger values of r.) The radiation causes the amplitude of the A-mode to 
decrease very slowly — by about 4% over 80 oscillations. We have done similar simulations 
at ,^ = 11 and ^ = 12 also, where we find ujq = 0.80 m and t^o = 0.98 m respectively. In 
these simulations, the oscillating soliton emits radiation with u = 2uJo, and the amplitude 
of the radiation and the rate of decay of the fundamental oscillation are larger than in 
Figure 5. The values of Uq for ^ = 10.4, 11, and 12 which we have found numerically are 
in good agreement with the relationship ( [4.16| ). This numerical evidence suggests that we 
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FIG. 5. In the left panels, we show p at r = 0.608/m as a function of time. It oscillates with 
period 16.69/m, and the amplitude of the oscillation is decreasing very slowly. In the right panels, 
we show p at r = lO/rn, to display the outgoing travelling waves shed by the oscillating soliton. 
These waves have three times the frequency of the fundamental oscillation seen at r = 0.608/m. 
Note that the amplitude of the outgoing waves is so small that they are invisible in the plots of 
p{r) on the preceding page. We conclude that for ^ = 10.4 the soliton has an almost stable mode 
of oscillation with frequency uq = 0.374 m — the A-mode — which slowly radiates waves with 
frequency Suq. 
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are justified in using the Hamiltonian (|4.13| ) to describe the long-hved normahzable A-mode 



with ujq < m and its couphng to the continuum. In the next chapter we will quantize this 
Hamiltonian and use it to describe the excitation of the A-mode by single IV-boson quanta. 

Finally we note that in principle it is possible to destroy a soliton with a minimum energy 
pulse, i.e. one whose energy is just above AE, and for ^ close to ^* this energy is small. To 
find the form of this pulse we could time reverse a solution which starts at the sphaleron and 
is given a gentle push towards the soliton. For ^ close to ^* so that the A-mode has a small 
frequency, the configuration takes a very long time to settle down to the soliton and in the 
process emits a very long train of low amplitude outgoing waves. Although the time-reversed 
solution consisting of a very long train of incoming low amplitude waves being absorbed by 
the soliton would eventually go over the sphaleron barrier and result in soliton decay, it 
would be rather difficult to set up initial conditions which produce this complicated, finely 
tuned, incoming configuration. Thus, the minimum energy soliton destroying pulses are not 
easy to build although we saw in Section II that with some extra energy, for ^ near C,*, the 
soliton is easily killed. 

V. QUANTUM PROCESSES IN THE FIXED AE LIMIT 

In the previous section we saw that for ^ close to ^* it is possible to identify a low 
frequency vibration of the soliton, the A-mode, with frequency ujq much less than m. If 
enough energy is transferred to this mode the soliton will decay. In this section we discuss 
the quantum mechanics of this mode. In this quantum setting the soliton can decay by 
barrier penetration as well as by being kicked over the barrier by a single VT-boson. We 
will see that if we work in a limit where AE is held fixed as we take g to zero, then we can 
reliably estimate the leading terms in both the tunnelling and induced decay rates. 

The Hamiltonian for just the A-mode coming from ([4.1 3|) is given by 
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where we have restored the g dependence. Note that ujq, b and all the terms in the ellipses 
depend on C, and m but not on g. We have changed the sign of A for later convenience. As 
^ goes to ujq goes to zero but the other terms are presumed not to change much. The 
classical soliton is at A = while the sphaleron is at \ = ujq/}) from which we have 

The fixed AE limit has g going to zero with ^ taken to ^* in such a way that ( ^.2|) is fixed. 
Since 6(^,m) does not vary much as ^ goes to ,^*, we see that in this limit ujq ~ (7^^^. Using 
( pl2D and ( [4.16| ), we see that g'^/^E ^ — ^*)3/2 go that in order to take the fixed limit 
we take g to zero with — ~ g^^^- (The reader who is concerned that the coefficient of 
A^ in uJqI g"^, goes to infinity in the fixed /S.E limit should note that because of the g^ 

in front of the in (5J) the frequency of oscillation is t^o-) 

When taking the fixed AE limit, it proves convenient to rescale according to 

\' = \uj^/g^\g-^l^ 

v' = pg/^o ^pg'^'^ 

h' = hg/ujlr^hg' . (5.3) 

Writing the Hamiltonian (|5.1| ) in terms of the new variables and then dropping the primes 
we obtain 

2 

= ^o' y + V{\) , (5.4) 



where 

V{X) = lx'-lx' + ... . (5.5) 

After rescaling, the sphaleron is at A = 1 /6 and the barrier height is given by 

AE = 1/(66^) . (5.6) 

Quartic and higher terms in V{X) are all suppressed by powers of g/too ~ g'^^^. Note that 
ujq now plays the role of h in the Hamiltonian (|5.4| ). As g goes to zero in the fixed AE limit, 

36 



ujq goes to zero like g^^^ and a semi-classical (WKB) treatment is appropriate in order to 
compute the leading sma\\-g behavior of the soliton destruction cross-section. 

In the fixed AE limit, the ground state of the quantum soliton has the A degree of 
freedom in a wave function V'o(A) which is described approximately by a harmonic oscillator 
ground state wave function: 

/ 1 N 1/4 / \ 2 \ 

There are three relevant scales in A, which differ in their (7-dependence. First, the width of 



the ground state wave function y (ipolX'^lipo) goes like y/uJo ~ g^^^. The second scale, which 
goes like g^, is the distance in A between the sphaleron at A = 1/6 and the minimum at 
A = 0. Note also that (|5.7| ) is a good approximation to tpQ for A such that the cubic term 
in V^(A) can be neglected relative to the quadratic term, namely for |A| ^ 1/6. Finally, 
note that the quartic and higher terms in V{X) can be neglected for A less than of order 
uJo/g ~ g""^^^, the third scale. Hence, as g is taken to zero with AE fixed, truncating the 
potential at cubic order becomes valid for larger and larger A. 

The soliton will decay if the A degree of freedom tunnels under the barrier given by the 
potential ^(A) shown in Figure 6. The rate is of the form 

r = (5.8) 

where the factor B is 

B.^r,,^W^^-W-'-l^^ (6.9) 

OJq Jo ^ 5 UqO^ 5 UJq 

We are able to neglect the width of the wave function ( |5.7|) in this calculation because as g 
goes to zero it is small compared to the change in A during the tunnelling process. Since in 
the fixed AE limit ujq ~ g^^^ we see that the tunnelling rate goes as exp(— constant/(7^/^). 
For the approximation to be reliable we require that B be much greater than one. This in 
turn requires that g be small. 

We can compare this calculation with that of Rubakov, Stern and Tinyakov [|1^ who 
numerically calculated the action of the Euclidean space solution which tunnels under the 
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FIG. 6. The potential V{X) for real A. For later use, the energies £"0 and E are also shown. 
■00 has three turning points, and A = Aq is the left-most of the three. tpE has one turning point at 
X = Xe. 



barrier. They used the equations of motion of the full 3 + 1 dimensional theory with the 
restriction to the spherical ansatz. At ^ = 12 we have AE = 1.2m/ g'^, ujq = 0.98m giving 
g'^B = 4.4 which is to be compared with what we read off Figure 2 of Ref. |]13|, namely 
g'^B = 4± 1. This agreement again supports the view that the A mode is the relevant degree 
of freedom for discussing soliton decay for C, near 

We now turn to induced soliton decay. Our picture is that the soliton will decay if the 
A-mode is excited to a state with energy above AE. The A-mode couples to the continuum 
modes which can bring energy from afar to the soliton. The free quantum Hamiltonian 
for the is 



dojoj 



2 2 I 

9^ 



ala^ + 1/2 



(5.10a) 
(5.10b) 



where 
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The qui have been chosen to diagonahze the fluctuation matrix at the sohton. Therefore 
Hq^ describes non-interacting massive W^-bosons propagating in a fixed sohton background. 
For each value of uo there are actually an infinite number of different VT-boson quanta. 
For example there are the states with frequency uj and aU values of angular momentum 
relative to the soliton center. These extra labels are omitted throughout but their presence 
is understood. 

The A-mode couples to the continuum modes through cubic couplings of the form 

Hmt = W\^l dud{u)q^ + —\ I duj e{uj , uj') q^q^X (5.12) 
which appear in (|4.13|) . We have rescaled A according to (|5.3|) . The couplings ( ^.12|) arose 



upon expanding about the soliton. The functions d{uj) and e(u;, oj') are peaked at values of 
UJ corresponding to wavelengths of order the size of the soliton. They are also only peaked if 
the unspecified labels allow large overlap with the soliton. For example even with uj chosen 
so that [uj"^ — m^)^^/^ ~ soliton size, it is only the low partial waves which have d{uj) and 
e{(jj, uj') large. 



The first term in (|5.12|) allows for the absorption of a single VT-boson by the soliton. 
The W^-boson energy E is transferred to the A-mode. The second term in ( |5.12|) allows a 
single W^-boson to scatter inelastically off the soliton, transferring energy E to the A-mode. 
We now calculate the rate for the absorption process; the calculation for the scattering 
process is similar. (The coefficients of the A and A^ operators have different (^-dependence, 
but this will not affect the leading ^f-dependence of the cross-section for either process.) 
Assuming that the soliton starts in its ground state, in order for the soliton to decay we 
require E + ujq/2 > AE. Since ujq <^ AE we can approximate this as E > AE. In the fixed 
AE limit we are free to choose AE to be a constant times m where the constant is of order 
unity. (Recall that m is held fixed throughout this paper.) Now the soliton size is roughly 
2/{my/^) and in the fixed AE limit ^ goes to ^* = 10.35. Thus the W^-boson wavelength 
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and the soliton size can be comparable. There is no length scale mismatch and d{E) need 
not be small. 

Using Fermi's Golden Rule we now calculate the cross section for iy + soliton anything 
with no soliton. Let |k) be a single W^-boson state with energy E, normalized to unit particle 
flux. Now 

(0|i/int|k) = 4 / dcodico) {0\q^\k)=g^d{k) , (5.13) 



where we have deflned (i(k) so that it is independent of g (see ( |5.11| )). The A- mode starts in 



the state iPq{X) with energy ~ C(Jo/2 which again we neglect relative to Ai?. The interaction 
( p.l3| ) can cause a transition to a state iPeW in which the A-mode has energy E. Since the 
width of -ipo is ~ g^^^ <^ 1, it is tempting to try approximating the states with E > AE as 
plane waves 

^^(^) ~ 1/2 Ov^AM) • (5.14) 

The cross section for a transition from ipQ to ipE is 



2 



O-dcstruction = A/" ^ I{E) , (5.15) 

where A/" is a ^f-independent constant and where I{E) is the integral 

X{E) = J dXM^)>^'M^) ■ (5.16) 
If we take ipo and ipE as in ( p.7|) and (|5.14|) respectively, X(i?) is easily evaluated, yielding 



J(E) ~ exp (-E/^o) , (5.17) 

where we have dropped all prefactors. This result is in fact incorrect .0 While it is true that 
( p.TD and ( |5.14] ) yield a good approximation to the integrand where the integrand is biggest. 



the result ( 5.17|) is exponentially smaller than the integrand. This raises the possibility 



^We are grateful to D. T. Son for noticing this, and for pointing us toward the correct answer. 



40 



that corrections to the wave functions neglected to this point may change (|5.17|) . We must, 



therefore, use WKB wave functions which take into account the quadratic and cubic terms 
in the potential V{X). As (7 — in the fixed AE limit, ujq and using semi-classical wave 
functions becomes a better and better approximation. We show below that for E = AE the 
leading dependence of the of as (? ^ in the fixed AE limit is in fact that of ( ^.17[ ) 
with the coefficient of AE/coq being (18 — 4-\/3)/5 instead of 1. Thus, we will find that even 
though the soliton destruction process does not involve tunnelling, the correct cross-section 
is exponentially small as cuq ~ g^^^ goes to zero. The reader who is not interested in the 
details of the evaluation of I{E) can safely skip to equation (|5.29|) . 
We now wish to evaluate the leading semi-classical dependence of 

IiEo,E) = I d\i,E,>?i^E (5.18) 

in the fixed AE limit where E > AE and AE > Eq > and where i/je and ipEo are WKB 
wave functions for the Hamiltonian ( |5.4| ). See Figure 6. The reader may be concerned that 
( p.l8| ) is infinite. (Both wave functions are real, and for large positive A the integrand ( |5.18| ) 



has a non-oscillatory piece which grows like A^A~'^/^.) However, when the relevant limits are 
taken correctly, the answer we seek is in fact finite. Recall that our problem reduces to that 
of the A mode in a cubic potential only for |A| < coo/g ~ g "^^^- Therefore, we should do the 
A integration from A = —A to A = +A, where A is real and positive and where we take A to 
infinity more slowly than g^'^^^ as g goes to zero. The result of such an evaluation would go 
like A'^/^ exp(— constant/ciJo)- Because we do not take A to infinity before taking g to zero, 
the prefactor does not make the result infinite. 

The evaluation of matrix elements of operators between semi-classical states has been 



treated by Landau |[T^, and although his final answer does not apply to our problem, we 
follow his method to its penultimate step. Landau's method yields only the leading [i.e. 
exponential) dependence of such matrix elements, and says nothing about the prefactors. 
Thus, using Landau's method yields the leading small-^f dependence of ( p. 18] ) irrespective 
of whether the prefactors make the integral infinite. In the calculation which follows, it 
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nevertheless proves convenient to multiply the integrand in ( ^.18|) by exp(— JA^/cuq) with J 
a constant. This does in fact render the integral finite, but it may also modify the exponential 
dependence of the result. Therefore, after the (7-^0 limit has been taken we must take the 
J ^ limit. Landau's method [0] applied to our problem yields 



I{Eo, E) ~ Im / d\ ^ 

y [{V{\)-E,){V{\)-E)f" 



X exp 



(^j^ dx^2 {y{x) - Eo) - dx ^2 - E) - JA^ ^ | 



. (5.19) 



In this expression, A is treated as complex and it is understood that the contour has been 
deformed into the upper half plane. This is done both in order to avoid the turning points 
on the real axis shown in Figure 7, and because in deriving (|5.19|) Landau uses expressions 



for WKB wave functions which are valid only in the upper half plane and not on the real 
axis. The first square root in the exponent in ( |5.19|) is taken to be positive on the real axis 



for A < Aq and the second is taken to be positive on the real axis for A < A^;. 

The equation V{x) — E = has three roots. One is at A^;, on the negative real axis, and 
the other two, at Abp and A^p, have nonzero imaginary parts. (For E AE, Abp goes to the 
real axis at Agph = 1/^-) evaluating ( |5.19| ) we must keep in mind that at A = Abp in the 



upper half plane, the integrand has a branch point. This singularity will play an important 
role in our analysis. (Unlike in the example treated explicitly by Landau, it does not arise 
from a singularity in V{X).) The branch cut from Abp must not cross the real axis, and it is 
convenient to take it to run upward vertically. The integrand in (|5.19|) is a function which 



is analytic in the upper half plane except at Abp and along the associated cut. To evaluate 
the integral, we are free to push the contour upward away from the real axis as long as we 
ensure that it does not touch the branch point Abp or cross the branch cut. 

We now evaluate the leading exponential dependence of (|5.19|) .P| To this end, we drop 
the prefactors in ( ^.19| ). We write the integral as 



^The analysis described below and the result ( |5.29| ) were provided by A. V. Matytsin. 
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'E '"0 

FIG. 7. This figure is a sketch showing the important points in the complex A-plane discussed 

in the text. The branch points are marked with dots, and the branch cuts are shaded. The contour 
in (5.19), just above the real axis, and the deformed contour we use to evaluate the integral are 
both shown. 



dX exp — \x{X)+iY{X) 



(5.20) 



where X and Y are real and where 
X + iY = 



dx ^2 {V{x) - Eo) - f dx ^2 {V{x) - E) - JX^ . 



(5.21) 



It is easy to check that for J = the integrand in ( |5.20|) has no saddle points at finite A. 
However, making J nonzero introduces a saddle point at large |A| which moves off to infinity 
as J is taken to zero and it is convenient for us to evaluate the integral with nonzero J and 
then take the J — limit. 

We now describe the behavior of X at large |A|. Write A = Aexpz6'. We have chosen 
the branch cut to run vertically and so it is at = 7r/2 for large A. To the right of the cut, 
that is for 7r/2>0>O, asA goes to infinity 



X ~ -A5/2gjj^(5^/2) - JA2cos(2^) 



(5.22) 
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and the J term is subleading. X goes to +00 at large A for 7r/2 > 9 > 27r/5 and goes to 
—00 for 271 / 5 > 9 > 0. The descent to —00 is most rapid for ^ = vr/S. To the left of the 
cut, that is for n > 9 > 7^/2, a.s A goes to infinity 

X ~ X* + A~i/2 sin(^/2) - JA^ cos(2^) , (5.23) 

where X* is a constant independent of J, A, and 9. (For J = 0, as A goes to infinity for 
TT > 9 > 7r/2, X X* and Y 0.) For nonzero J, there is a saddle point at finite A. For 
small J, this saddle point is at ^ ~ Svr/S and A ~ J^^/^. Thus, as J the saddle point 
recedes to infinity as promised, and JA^ at the saddle point goes to zero. Therefore, in the 
J — s> limit X at the saddle point goes to the value X* . 

We now deform the contour as sketched in Figure 7. For nonzero J, the saddle point is 
at finite A and we choose the contour to follow the path of steepest descent from this saddle 
point. To the left of the saddle point, the steepest descent path curves toward the real axis, 
and then approaches the real axis asymptotically. As we discuss below, X(Abp) is greater 
than X* . Therefore, to the right of the saddle point, the path of steepest descent from the 
saddle point cannot get around the branch point and necessarily runs into the branch cut. 
After reaching the cut, the next section of the path ascends as it traverses (II), following the 
cut inward toward the origin, until it reaches the region of the branch point Abp. Along (II), 
X ascends monotonically from X* to X(Abp). Y is not constant. Then, to the right of the 
cut, the contour follows the path of steepest descent (III) toward infinity along 9 = tt/S. 

There are two contributions to the integral ( |5.2(j| ). First, the saddle point makes a 
contribution which goes like exp(X*/ci;o)- (Note that we take the g ^ limit and then 
take the J —>■ limit.) The second contribution arises because the path must ascend from 
the saddle point at infinity as it traverses (II) in order to get around the branch point, 
before then descending along (III) to the right. Therefore, the integral ( |5.20|) receives a 
contribution from the region of the branch point which goes like exp{X {X\yp) / ujq) . In sum, 
therefore, the integral ( p. 191) goes like 



I{Eo, E) ~ exp(X7a;o) + exp(X(Abp)M) , (5.24) 
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where we have dropped the prefactors, about which Landau's method says nothing. At this 
point, we can take the Eq ^ and E —>■ AE hmits simply by setting Eq = and E = AE. 
Prior to this point in the calculation, taking these limits would require careful treatment of 
branch points. Henceforth we set Eq = and compute T{E) = X(0, E). 

It only remains to evaluate the relative size of X(Abp) and X*. Both X(Abp) and X* 
depend on E. After some calculation one finds that for E = AE 

X* = -uoB = -^AE , (5.25) 

where B is the tunnelling amplitude computed in ( |5.9|) , and 

X(A.,) = -v/2f' dAv/AV2-.AV3 = -^(3-A) . Jl^E - ^) . (5.26) 

SO X(Abp) is the larger {i.e. least negative) of the two at = AE. At large E, both 
X(Abp) and X* decrease like —E^^^. For E > AE, the integrals in ( ^.21|) must be evaluated 
numerically. We find that both X(Abp) and X* decrease monotonically with increasing 
energy, and X(Abp) is always greater than X*. Consequently, the integral is dominated by 
the region of the branch point for all energies E > AE. That is, 

J(E) ~ exp(X(Abp)M) (5.27) 

and 

f destruction ~ (P(k)exp(2X(Abp)M) , (5.28) 

where we have dropped all prefactors except d. Thus, although the integrand has a saddle 
point (at infinity), the integral is not dominated by that saddle point. This occurs because 
the path of steepest descent from the saddle point necessarily runs into the branch cut. 
Equivalently, the presence of the branch cut prevents the actual contour of integration from 
being deformed into a path of steepest descent through the saddle point. Although the path 
can be deformed to pass through the saddle, it must ascend from the saddle to the region 
of the branch point. (Note that although X(Abp) > X* for all energies E > AE, X(Abp) 
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is greater than Buq, and the rate for induced sohton decay is greater than the tunnelhng 
rate, only for E within a range of energies which we determine numerically to be AE < E 
< 1.7AAE.) 

Because X(Abp) decreases monotonically with increasing E, the cross section ( p.28| ) for 



the soliton to be destroyed by a single VT-boson is least suppressed by T{E) at threshold. 
For E = AE the soliton destruction cross section goes like 



CTdestmction ~ d (k) CXp I AE / UJq I 



(5.29) 



as — > in the fixed AE limit. 

We expect d{E) and accordingly (i(k) to be appreciable when E ~ AE so long as 
AE is comparable to the inverse soliton size, which is of order the inverse VT-mass. Un- 
der these conditions, there will be no length scale mismatch and d{E) will not depend 
sensitively on E for E ~ AE, so o-dcstruction will be maximized for E = AE. Thus the 
maximum rate for soliton decay induced by collision with a single ly-boson is proportional 
to exp(— (36/5 — 8^/3 / 5) AE / luq) . This is to be compared with the tunnelling rate in the 
same limit which is proportional to exp(— (36/5)Ai?/u;o). Both go to zero as g goes to zero 
like exp(— constant/^f^/^), but the ratio of the tunnelling rate to the induced decay rate is 
exponentially small. 

We have computed the cross section for a single PF-boson to be absorbed by the soliton 
and to excite the A-mode to a continuum state above the barrier, which in our picture results 
in soliton decay. The cross section for a tV-boson to destroy the soliton by scattering off the 
soliton and transferring energy E to the A-mode can be calculated using the second term 
in ( ^.12|) . The calculation is similar to the one we have done and the result has the same 
exponential factor as in ( [5.29| ) but would have a different prefactor. Because the exponent 



in (|5.29| ) includes cOq'^ ~ g~^^^, these cross sections go to zero faster than any power of g 
as g goes to zero in the fixed AE limit. Note that this suppression arises even though the 
process does not involve tunnelling and even though there is no length scale mismatch. It 
arises as a consequence of the limit in which we have done the computation, because in that 
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limit destroying the soliton reduces to exciting a single degree of freedom to an energy level 
infinitely many (~ AE/ujq) levels above its ground state. Thus, taking g at fixed AE 
makes the computation tractable but makes the induced decay rate exponentially small, 
albeit larger than the tunnelling rate. 



We have described classical and quantum processes in which electroweak solitons are 
destroyed. In this section, we argue that if we couple a quantized chiral fermion to the gauge 
and Higgs fields considered in this paper, then soliton destruction implies nonconservation 
of fermion number. The argument we present treats the gauge and Higgs fields as classical 
backgrounds. In particular, we ask how many fermions are produced in a background given 
by a solution to the classical equations of motion in which a soliton is destroyed. We expect 
that our conclusions will also be valid for soliton destruction induced by a single W^-boson. 

We introduce a quantized fermion field \E', and as in the standard electroweak theory but 
neglecting the U{1) interaction, we couple only the left-handed component of the fermion to 
the non-Abelian gauge field. We add the usual Yukawa coupling between the fermion and 
the Higgs field to give the fermion a gauge invariant mass. The Lagrangian for the fermion 
is 



where D^ = d^- iA^Pl, = i(l - 75) and = |(1 + 75). The Higgs field $ of Q is 



VI. FERMION NUMBER VIOLATION 




(6.1) 



given by $ = {v/^/2)U. For simplicity, both the up and the down components of \E' have 



the same mass rrif. The gauge invariant normal ordered fermion current 



Jf" = : ^7^*^ : 



(6.2) 



is not conserved, that is. 



1 



(6.3) 



327r2 
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We consider backgrounds given by solutions of the kind found numerically in Section II. 
After the soliton has been destroyed the solution dissipates. By dissipation we mean that 
at late times the energy density approaches zero uniformly throughout space. This means 
that at late times the solutions are well described by solutions to the linearized equations of 
motion 

(d^d" + m^) Al"" = (6.4) 



in unitary gauge. It is tempting to try to integrate ( |6.3| ) and relate the fermion number 
violation to the topological charge 

Q = ^2 fd'^ '"'''^ Tr . (6.5) 

First, note that because the region of space-time in which Ff^,y 7^ is not bounded, there 
is no reason to expect Q to be an integer. Furthermore, it is shown in Ref. |jl9| that for 
a background which satisfies ( |6.4| ) at early and/or late times the integral in ( |6.5| ) is not 
absolutely convergent and Q cannot sensibly be defined. 

In a background given by a solution to the equations of motion which dissipates at early 
and late times, the number of fermions produced is known to be given by the change in 



Higgs winding number • In this paper, the Higgs mass is infinite so the Higgs winding 

number can never change. For solutions with no solitons in the initial or final states, the 
arguments of Ref. apply, and no fermions are produced. However, if there is a soliton 



in the initial or final state the assumption of Ref. [|T9[ that the solution dissipates is not 
satisfied. In this section, we show that in a background given by a solution in which one 
soliton is destroyed, one net anti- fermion is produced if the fermion is light {rrifL <^ 1 
where L is the size of the soliton) and no fermions are produced if the fermion is heavy 
[rrifL ^1). In the rufL ^ 1 case, however, there is still a violation of fermion number in 
the sense that the soliton carries heavy fermion number whereas the dissipated configuration 
after the soliton is destroyed does not. 

We now review some known facts about fermion charges which a background field con- 
figuration can carry |pll,p|,p2| . Consider some localized time independent field configuration 
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yl^(x) in the unitary gauge U = 1. Imagine adiabatically interpolating from the trivial 
background A^{x.) = to the A^{x.) of interest and following the adiabatic evolution of the 
fermion state which at the beginning of the interpolation is the fermion vacuum. At the 
beginning of the interpolation, the state has all negative energy levels filled with the mode 
functions determined by the single particle Dirac Hamiltonian 

i^ferm(t) = 7° [-^7^ A + nif] - AqPl (6.6) 

with = 0. The change, from the beginning to the end of the interpolation, of the expec- 
tation of the fermion charge operator, / rf^x J*^ with J° given by (|6.2|) , has been calculated 



| 21| . The result is the Chern-Simons number of the field A^ 



3 

AiAjAk + — iFijAk 



(6.7) 



Since at the beginning of the interpolation the fermion charge is zero, A^cs is the charge of 
the state arrived at by adiabatically following the initial vacuum state. We will see that the 
fermion state reached by this adiabatic process is not necessarily the lowest energy fermion 
state in the background A^(x). 

Consider the case when the final configuration has the special form Aj(x) = iUldiUi, 
Aq = where Ui is a winding number one map, say of the form (|2.9|) , with a characteristic 
size L. In this case, A/'csf^] is the winding of Ui, that is Ncs[iUidiUi] = 1. Thus the state 
arrived at at the end of the interpolation has fermion charge one. We now examine what 
this state is. At all times during the interpolation we can define an instantaneous single 
particle Dirac Hamiltonian by ( |6.6|) and we can therefore discuss how the spectrum of the 
instantaneous Hamiltonian varies during the interpolation. At the beginning we have the 
free massive Dirac Hamiltonian which has a gap between —mf and rnj. If rrif ^ 1/L, 
then the spectrum is not perturbed much by the gauge field and in particular no energy 
levels cross zero during the interpolation. In this case, throughout the interpolation the 
fermion state is the lowest energy fermion state in the presence of the bosonic background. 
However if rrifL <^ 1, it has been shown that one level crosses zero from below during 
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the interpolation. This means that for rufL <^ 1, the state reached at the end of the 
interpolation is the lowest energy fermion state in the presence of Ai = iUldiUi plus a single 
fermion. Thus the charge of the lowest energy fermion state in the presence of Ai = iUldiUi 
is zero for m/L ^ 1. However for nifL ^ 1, the charge of the lowest energy fermion state 
in this background is one. This ends our little review. 

In the case at hand the electroweak soliton ^4^°' is only approximately of the form iUldiUi 
and the Chern-Simons number of the soliton configuration is not an integer. In fact a general 
background A^{:x) will not carry an integer- valued fermion charge and this charge can be 
viewed as a consequence of the polarization of the vacuum by the background. Nonetheless 
we will argue that when the soliton is destroyed an integer number of fermions is produced. 

Consider a background given by a solution in which a soliton is destroyed. At t = — Tq ^ 
the solution consists of a soliton at rest and an incoming pulse and at t = Tq ^ there 
is only outgoing radiation. We wish to avoid evaluating the fermion charge aXt = iTo. To 
this end we introduce a background configuration A^{^,t) for —T < t < T with T > Tq 
which agrees with the solution A^(x, t) for — Tq < t < Tq. At t = — T we choose the 
background A^(x, — T) = iUldiUi where Ui is a winding number one map which produces a 
configuration which is close to the soliton, that is, A^°^ ~ iUldiUi. In particular the length 
scale over which Ui varies, L, is determined by the size of the soliton. Thus the interpolation, 
running backward from — Tq, turns off the incoming pulse and distorts the soliton until it 
is of the form iUldiUi. Now at t = Tq the solution is just the outgoing remnants of the 
soliton and the initial pulse, and at t = T we chose A^(x, T) = 0. Thus between Tq and 
T the interpolation turns off the outgoing pulse bringing the background to its vacuum 
configuration. 

The interpolation v4^(x, t) begins at t = — T with a configuration of the form Aj(x, — T) = 
iUldiUi and ends at y4^(x, T) = 0. The scale of Ui is L. Consider a fermion field coupled 
to this background with rrifL ^ 1, the heavy fermion case. The initial gauge field configu- 
ration has heavy fermion number one. Consider the instantaneous Dirac Hamiltonian ( |S.6[ ). 
Throughout the interpolation the gauge field is small compared to m/ and we conclude 
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that no levels cross zero throughout the interpolation. Thus even without a detailed field 
theory description of fermion production we conclude that the fermion state we reach at the 
end of the interpolation has no extra heavy fermions. The final configuration is v4^(x, T) = 
so the fermion number of this gauge field configuration is zero. The heavy fermion number 
of the gauge field background has changed from one to zero and no heavy fermions have 
been produced and so we see an anomalous violation of heavy fermion number. 

We now turn to the light fermion case, nifL <^ 1. At t = — T the configuration 
ylj(x, — T) = iUidiUi has light fermion number zero and we begin in the lowest energy 
fermion state in this background, which has all negative energy levels filled and all positive 
energy levels empty. Both the light and the heavy fermion number currents have the same 
anomalous divergence ( |6.3| ) so the difference between light and heavy fermion numbers is 
strictly conserved. We conclude that at t = T the state we arrive at must have light fermion 
number minus one. The final configuration y4^(x, T) = has light fermion number zero. 
Thus we see that the fermion state at t = T has one more light anti-fermion than light 
fermion. Thus in the background between t = —T and t = T no heavy fermions are 
produced and one net light anti-fermion is produced. 

We have discussed fermion production in a background going from A,(x, -T) = iUldiUi 
to Aj(x, T) = whereas we are really interested in fermion production only in the presence 
of the solution y4^(x, t) for — Tq < t < Tq. Therefore we need to argue that no fermions 
are produced (light or heavy) between — T and —Tq and between Tq and T. By making Tq 
arbitrarily large, we can make the amplitude of the incident and outgoing pulses arbitrarily 
small, and so make their effect on the Dirac Hamiltonian arbitrarily small. This ensures that 
no fermions are produced during the interpolation between Tq and T. It also ensures that, 
working backwards in time from —Tq, we can interpolate to a configuration with = A^°^, 
removing the incident pulse, without producing any fermions. 

It only remains to consider the interpolation between v4j(x, — T) = iUldiUi and A^ = 
A^°^. We can choose Ui to be the winding number one map which characterizes the skyrmion 
with the same e and v as the soliton of interest. We can then choose the interpolating 
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configurations to be a sequence of solitons witli fixed e and v witli g clianging from tlie value 
of interest to zero. Tlie beliavior of x during such an interpolation is depicted in Figure 2. 
Note that in taking g to zero with fixed e and v, ^ goes to infinity and m goes to zero in such a 
way that the soliton size stays fixed. Note also that while we are choosing the configurations 
during the interpolation to be solitons with differing values of g, the coupling between the 
gauge field and the fermions is held fixed. For frifL ^ 1, throughout the interpolation the 
fermion spectrum is not perturbed much by the gauge field, and so no levels cross zero. 
In our little review we learned that the configuration iUldiUi has fermion number one if 
nifL ^ 1 and fermion number zero if TUfL < 1. Since Ncs[iUldiUi] = 1 this means that 
as we reduce rrif from ruf ^ 1/ L to rrif ^ 1/L one (net) level must cross zero from below. 
Accordingly, for one or more values of rrif, of order 1/L, the Dirac Hamiltonian (|6.6| ) has 
a zero energy bound state in the skyrmion background and furthermore there is a nonzero 
value oi rrif below which there are no zero energy bound states. Now consider interpolating 
from the skyrmion configuration at t = — T to the soliton configuration. Define m*f{t) to 
be the largest value of rrif such that for all < m*f{t) the Dirac Hamiltonian ( |6.6| ) in the 
background yl^(x, t) does not have a zero energy bound state. From our understanding of 
the skyrmion background, we know that m*f{—T) is of order 1/L. By making C, arbitrarily 
large, the difference between the skyrmion and soliton configurations can be made arbitrarily 
small, and accordingly the change in the spectrum of the Dirac Hamiltonian during the 
interpolation can be made arbitrarily small. Therefore, for large enough ^, throughout the 
interpolation from the skyrmion to the soliton m*f{t) remains nonzero. We now assume that 
this is in fact the case for all > .^*. We feel that this is a reasonable assumption since, 
as Figure 2 shows, the soliton configuration is quite similar to the corresponding skyrmion 
configuration even for ^ near Making this assumption, we conclude that for m^L ^ 1, 
as for rrifL ^ 1, no levels cross zero during the interpolation between v4j(x, — T) = iUldiUi 
and = A^°^. Hence, no fermions (light or heavy) are produced during the interpolation 
between — T and — Tq. Therefore, in the background between —To and Tq, which is a classical 
solution in which a soliton is destroyed, no heavy fermions are produced and one net light 
anti-fermion is produced. 
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Suppose we are only interested light fermion production. We can view the heavy fermion 
as a device introduced only for the purpose of making an argument. Because we have not 
included the back reaction of the fermions, heavy or light, on the bosonic background, any 
conclusions we reach about the light fermion are in fact independent of whether there is or is 
not a heavy fermion in the theory. Therefore, in any process in which a soliton is destroyed, 
one net anti-fermion from each light SU{2)l doublet is anomalously produced. 

VII. CONCLUDING REMARKS 

We have described a theory which agrees with the standard electroweak model at 
presently accessible energies but which includes a metastable soliton with mass of order 
several TeV. This Higgs sector soliton may have a dual description as a bound state particle 
made of more fundamental constituents or it may be that the Higgs sector is fundamental 
and when quantum effects are taken into account, a metastable soliton is found. In any 
event, given the soliton, under certain circumstances we can reliably estimate the rate for 
collision induced decays. The parameters of the theory can be chosen so that the soliton 
configuration is close to the sphaleron configuration, which means that using the soliton as 
an initial particle makes it easy to find sphaleron crossing processes. Indeed, we have found 
classical solutions in which the soliton is destroyed where the incoming pulse corresponds 
to a quantum coherent state with ~ 1^-bosons. The rate for such processes is not 
exponentially suppressed as g goes to zero. Furthermore in the limit g goes to zero with 
— Mgph — Mgoi fixed we can reliably estimate the rate for a two particle scattering 
process in which a single incident VF-boson kicks the soliton over the barrier causing it to 
decay. We have argued that in all processes in which the soliton disappears fermion num- 
ber is violated. This model may be relevant only as a theoretical foil, as a demonstration 
that fermion number violating high energy scattering processes can be very different than in 
the standard model. However if no hght Higgs boson is discovered, it is even possible that 
Nature may be described by such a model. 
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APPENDIX A: 

In Section II, we presented the equations of motion in the spherical ansatz in the unitary 
gauge. In this appendix, we begin by presenting the action and equations of motion in the 
spherical ansatz without fixing a gauge. All the numerical solutions presented in this paper 
were obtained by solving both the unitary gauge equations of motion and the Aq = gauge 
equations of motion using different numerical schemes. We sketch both methods in this 
appendix. 

The spherical ansatz is given by expressing the gauge field and the Higgs field $ in 
terms of six real functions oq , ai , a , 7 , /i and u of r and t. is given in ( |2.3| ) and $ is 
given by 

<l>(x,t) = ^ [/i(r,t) + iz/(r,t)a- x] . (Al) 

For A^ and $ to be regular at the origin, we require that ao, a, ai — a/r, 7/r and z/ vanish 
as r 0. It is convenient to define the complex field 
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(A2) 



When the Higgs mass is set to infinity, |</)| is frozen at its vacuum expectation value ^/2m 
and 



m exp ir]{r, t) 



(A3) 



with 1] vanishing at the origin. Under a gauge transformation of the form exp[iQ{r, t)a ■ x/2] 
with i7(0, t) = 0, configurations in the spherical ansatz remain in the spherical ansatz and 
continue to satisfy the appropriate boundary conditions at the origin. Thus, the SU{2) 
gauge theory reduced to the spherical ansatz has a residual U{1) gauge invariance. 

In the spherical ansatz the action associated with the Lagrangian (|1.8| ) takes the form 



S 



Air 



oo 

dt I dr 



1 
2^ 



\X\ 



2m^ Re (ix^e'*") 



-4 9^r/ 



2r 



(A4) 



where 



X + ie 



2iri 



(A5) 



and the rest of the notation is as in Section II. The notation is chosen to manifest the U{1) 
gauge invariance present in the action ( [A4D . The complex scalar fields x ^i^nd have U{1) 
charges of 1 and 1/2 respectively, is the U{1) gauge field, f^u is the field strength, and 
is the covariant derivative. The indices are raised and lowered with the 1 + 1 dimensional 
metric ds"^ = dt^ — dr"^. 

The equations of motion in the spherical ansatz are 



1 



(A6a) 



i\x\ 



-m 



X 



1 / 

4e(^ + 



2iri 



1 \2 1 



(A6b) 
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Re(xV-) + ^Re(xV- 



(A6c) 



The same equations are obtained either by varying the action ( p..8| ) and then imposing the 
spherical ansatz or by varying the action (|A4|) . The unitary gauge action and the equations 
of motion of Section II are obtained by setting r] = 0. As a consequence of the U{1) gauge 
invariance, the five equations above are not independent, and in Section II we chose to 
discard (|A6c|) . 



In presenting the results of Section II, we found it useful to use the variables p, 9, and ai. 
The first is explicitly gauge invariant. The latter two are specific to the unitary gauge, and 
are equivalent to the gauge invariant variables {6 — 2r]) and (ai — 2diri). For an extensive 
discussion of gauge invariant variables in the spherical ansatz, see Ref. ||23|| . 

In the simulations of Section II, the boundary conditions at the large-r boundary of 
the lattice are not important, since we are interested in an ingoing pulse and its effects 
on the soliton and we stop the simulation before the remnants of the soliton reach the 
large-r boundary. In doing the simulation in Figure 5 of Section IV, we must have energy 
absorbing boundary conditions at the large r boundary of the simulation so that we can 
see the amplitude of the A mode oscillation decrease as it radiates energy away. At large 
radius, x satisfies [d^/dt"^ — d'^/dr^ — m^]x(r, t) = 0. We choose to impose x = —dx/dr 
at the large r boundary. This has the virtue that the energy flux at the boundary is never 
negative — the boundary can only absorb energy and cannot emit it. 

We now give brief descriptions of the methods we used to solve the unitary gauge and 
Aq = gauge equations. In Aq = gauge, we chose to write a lattice version of the action 
( |A^ ) and vary it, thus obtaining discretized versions of the equations of motion. The fields 
X and rj live at the sites of the lattice, while Oi lives on the spatial links. For a more detailed 
description of this technique as applied to the 1 + 1 dimensional Abelian Higgs model, see 
Refs. ||2^J2^ . In the Aq = gauge, we have the freedom to make a time-independent gauge 
transformation and we used this to set //(r, 0) = 0. However, rj does not remain zero at later 
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times. We specified initial conditions for x, x, ai and ?) at t = 0, and tlien cliose initial 
values for hi such that Gauss' law is satisfied. The equations of motion are second order in 
time derivatives of the fields 0,1 and rj. We used the fields at two successive time steps to 
determine the fields at the next time step. We used a lattice with 4000 sites and a lattice 
spacing dr = 0.004/m, and used a time step dt = dr/2. Because the equations of motion 
were obtained by varying a lattice action, there is a lattice version of Gauss' law which is 
satisfied during the evolution up to the precision allowed by computer arithmetic. Energy 
conservation was satisfied throughout the evolution to better than one part in a thousand. 
We verified that reducing dr and dt did not change any results quantitatively, and reduced 
the violation of energy conservation. We also checked that if we take the final configuration 
from a simulation like that of Figure 3, change the sign of all time derivatives, and run the 
simulation backwards in time, we obtain the initial configuration of the original simulation. 

With the unitary gauge equations of motion, we followed a different strategy. We first 
wrote them as first order equations for oq, ai, and the canonical momenta Ua^ and 11^. Of 
these seven equations, one is redundant. In doing the numerical evolution it was convenient 
to drop the equation for Hq^. We therefore chose initial conditions by specifying x, ai, oq, 
and n^, and then chose such that Gauss' law was satisfied. Rather than using a lattice 
action, we discretized the equations of motion themselves. We took care, however, to put 
X, and ao at the lattice sites while putting ai and Ua^ on the links, and to discretize 
the equations of motion in such a way that all quantities were accurate to one order in dr 
beyond the trivial one. We implemented the time evolution using a fourth order Runge- 
Kutta algorithm. We worked with a lattice spacing of dr = 0.02/(mv^), equal to 0.0058/m 
for ^ = 12, and a time spacing dt = 0.002/m. We verified that energy conservation and the 
equation for Ila^ were satisfied. For smaller dr and dt, no qualitative changes occurred and 
both checks were satisfied more accurately. 

APPENDIX B: CLASSICAL SOLUTIONS AND COHERENT STATES 

There is clearly a relation between classical solutions in Minkowski space and the tree 
approximation in quantum field theory in this section we derive a version of this relation 
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useful for our purpose. This appendix is self-contained and can be read independently of 
the rest of this paper. 

For simplicity, consider real scalar field theory with the action 



5 = -| d^xj- 9^09^0 --mV-V^(0)} (Bl) 
where V{(j)) is a polynomial containing cubic and higher terms. The classical field equation, 

+ m^) + V'{(t>cix)) = (B2) 
is equivalent to the integral equation 

= F^x) + F:^,ix) - 1 1 Df{x - y)V'[<i>,{y)Yy (B3) 



where 



1/1 tJ^Sj 



is the Feynman propagator and 



F,^ix)= dkf,^{k)e-''^ . (B5) 

out out 



In (|B5| ), = \\^W + rn?\ and dk is the Lorentz invariant measure 

dk = — ^ 4?^ d^k (B6) 

We do not yet require that (pc{x) is real, /in and /out are then arbitrary complex functions 
of k. 

Equation ( P3|) can be solved by iteration, and the solution expressed as a sum of terms 
described by Feynman diagrams. A typical diagram (for V{(f>) = ^(f)^) is shown in Figure 
Bl. Its contribution to (pdx) is 

yd'yd'z{-zfDF{x-y)DF{y-z){F:,,{z)}'FM . (B7) 

The sum is over all connected, tree diagrams. The combinatorial factors are determined by 
the symmetry of the diagrams in the usual way. 
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FIG. Bl 



We now ask which quantity in the quantum field theory with action ( pTl) is given in tree 
approximation by the same set of diagrams. The answer is 

(/out|0(a;) I /in) 



M^) = {<P{x))f.. 



tree 



(/out I /in) 



[B8) 



tree 



where 



/in ) are coherent states defined by 

out' 



1 



out 



/in ) = ^77 / /in (ki)dki ... /in (k„)dk„ ki . . . k„; 

n 9^- out 



out 



out 



(B9) 



where 



ki . . . k„; out ) are the usual asymptotic n-particle states. In terms of the quantum 



field (h(x) 



f 



in 

out 



. 



(BIO) 



Note that with our definitions of S and |/), both sides of ( P8|) are independent of g. The 



one-loop corrections to are of order g^. Note also the normalizations 



I = E ;^/dki . . .dk„ |ki . . .k„) (ki . . .k„| 



(Bll) 



and 



{mx)\0)=ge 



ikx 



fi 



out 

(/in /in ) = exp 



out 



out 



9^ J c 



m 
out 



dA; . 



(B12) 
(B13) 
(B14) 



It is important to note that the number of particles in the states |/) is of order l/g"^. 

There is certainly a large class of solutions ipdx) which behave at early and late times 
like incoming and outgoing wave solutions of the free field equation, with amplitudes which 
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go to zero at early and late times because of both the spread in 3-space and the dispersion 
due to m 7^ 0. For such solutions, with /j^ smooth enough, we have 



out 



-oo, (j)c{x) ~ J dA;/in(k)e-^^^ + J dfc /i*„(k)e^*^^ (B15a) 
^ +00, <j),{x) ^ Jdk Knt{k)e-"'^ + J dk f^Me'"' (B15b) 



where, using (B3) and the values of Df{x — y) as — > =Fc>o, 



h*M = /:ut(k) - ^ J d%e-'''yV'{Uy)) (B16a) 

^out(k) = /in(k) - t J d'ye^'^yV'^My)) ■ (B16b) 

With a given complex (pd^), we read off /in(k) from the positive frequency part as x^ —00, 

and /*^^(k) from the negative frequency part as x^ +00 to discover what (pdx) is the tree 
approximation to. If <f)c{x) is real, 

/ii,(k) = /i„(k) (B17) 

out out 



and equations ( [B16[ ) become equivalent relations which determine /out(k) in terms of /in(k) 



Of course we are concerned with given real solutions (pd^) and we simply read off both 
and /out from the asymptotic behavior; we do not need equations ( [B16D to relate them. 

Returning to the complex case, we derive an expression for the matrix element (/out I /in) 
in the tree approximation in terms of the corresponding 0c- The quantum perturbation 
theory gives 

(/out|/in)=explT(/o,t,/in) (B18) 

9^ 

with T expressed as the sum of connected Feynman diagrams. A typical tree diagram (again, 
for V{(f)) = h4>^) is shown in Figure B2. This contributes 



6 



J Dp{x - y)DF{y - z) {F,^{x)f F;,,(|/)Fi„(^)F:,J^)d^x d^y dS (B19) 



to T{fo^t, /in). Now consider small changes 5/in . It can be seen that 
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FIG. B2 



SnUt, /in) = J dk6fUk) |/i„(k) + z J d'ye^'y {d' + m') {4>{y))f 

+ J dk6Uk) {/:„t(k) + t J d'ye-^'y {d' + m') (B20) 



where 



(/out|0(a;) I /in) 



(B21) 



(B22) 



(/out I /in) 

In the tree approximation, using ([BHD, ( P^ ) and ( |B16| ), 

5Tt,ee(/out,/in) =/dfc {<5/:„, (k) ^out (k) + 5/in(k) /l^nlk) } • 

We can now verify that ( |B22| ) is satisfied by 

Ttree(/out, /in) = ^ / /o^nt (k) /lout (k) + ^ / fMhU^) 

-I I d^x l^Mx) {{d^ + m2)0e(a;)} + V^(0c(a;))] . (B23) 

Of course ( P23|) is i times the classical action plus boundary terms. The point of the 
preceding argument is to get the boundary terms right. The expression ( P23|) can also be 
derived from the stationary value of a functional integral expression for (/out | /in)- Note that 
the integral over space-time can be written as 



e 



d^x 



V(0,(x)) - h,{x)V'{Ux)) 



(B24) 



and so converges for typical solutions (pc{x) which fall off rapidly enough in both space and 
time. 

Now consider /jn which correspond to a rea/ solution 4>cix). Then (using ( |B18| )), 



out 
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(/out I /in) 



tree 



exp 



^/l/o..rd* + 5^/l/.nl=d*- 



where 6 given by ([B24|) is real. Using (|B14|), we see that 



7^ 



(B25) 



I (/out I /in) 



tree I 



1 . 



( /out I /out ) ( /in I /in ) 

One-loop corrections to T(/out, /in) are of order g"^. Thus if we consider states 



(B26) 



/ 



in 
out 



/i 



out 



/in / 



m 

out out 



1/2 



(B27) 



normahzed to unity, 



/J=e-^^/^\ /out) + 



(B28) 

where /out) = and \a\ is not of order exp(— c/(7^) as — 0. In fact (not proved here) 
to order gf^, i.e. one loop, 



I |2 -P 

a = e 



(B29) 



where P is the probability of the transition vacuum one pair of particles in the theory 
with action 



S' = ]- ( d^x {9^09^0 - m2(x)02} 



(B30) 



where 



m^{x) = m'^ + V"{(f)cix)) . 



(B31) 



Thus, to each real solution (pd^) to the Minkowski space classical equations of motion, 
there corresponds a not exponentially suppressed scattering process between coherent states 
(containing of order l/g"^ particles). 
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